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Abstract 

Waves propagating inwardly to the wave source are called antiwaves which have negative phase 
velocity. In this paper the phenomenon of negative phase velocity in oscillatory systems is studied 
on the basis of periodically paced complex Ginzbug-Laundau equation (CGLE). We figure out a 
clear physical picture on the negative phase velocity of these pacing induced waves. This picture 
tells us that the competition between the frequency to ou t of the pacing induced waves with the 
natural frequency Uq of the oscillatory medium is the key point responsible for the emergence of 
negative phase velocity and the corresponding antiwaves. u> out uJo > and \u> ou t\ < \u>o\ are the 
criterions for the waves with negative phase velocity. This criterion is general for one and high 
dimensional CGLE and for general oscillatory models. Our understanding of antiwaves predicts 
that no antispirals and waves with negative phase velocity can be observed in excitable media. 
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I. INTRODUCTION 



A kind of novel waves, waves propagating inwardly to the corresponding wave sources, 
were predicted theoretically in sixties of the last century [l|, and were experimentally dis- 
covered at the beginning of this century 2], [3]. Throughout this paper, we call this spe- 
cial kind of waves antiwaves (AW) for short, which have negative phase velocity, and call 
waves propagating forwardly from the wave source normal waves (NW). The phenomenon 
of negative phase velocity exists in both linear^ and nonlinear [b| systems. So far, lin- 
ear optical antiwaves have attracted great attention and have been introduced to applied 
sciences 101 ] . While nonlinear antiwaves have been much less. Till now, 

scientists have found only antispirals in nonlinear systems [ll. 12, 13, Q, [l^. In contrast 



to normal spiral waves, antispiral waves propagate inwardly to the spiral tips, which serve 
as the sources of the waves. In some recent works, several models have been suggested to 
understand antispirals 0, ll|. Gong and Christini have put forward an empirical criterion 
for the emergence of antispirals |l3]. Bar et. al. have explained antispirals from the changes 
of the sign of u (spiral frequency) and k (spiral wave number) with the model of CGLE 16]. 
While these studies illuminated some aspects of antispiral waves, the mechanism and the 
general physical picture of antiwaves are still not completely clear. With antispiral waves 
in mind, some questions emerge naturally: Are there general antiwaves exist (beside anti- 
spirals) which propagate toward the wave sources? If yes, is there any common mechanism 
underlying these antiwaves? 

This paper will deal with the above problems. The paper is organized as follows. In 
Section 2, we specify our nonlinear model to produce travelling nonlinear AWs by certain 
local periodic pacing. In Section 3, we analyze the mechanism of general nonlinear AWs, 
including antispirals and planar travelling antiwaves as their special cases. In Section 4, we 
specify the parameter regions for emergence of AWs in CGLE. The CGLE system serve as 
a simple example showing applications of our understanding on nonlinear AWs. Since our 
analysis does not depend on the detail dynamics of CGLE, we expect that the analysis in 
this paper is generally valid for oscillatory media. We conclude our results and give some 
remarks in Section 5. 
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II. MODEL AND NONLINEAR WAVES WITH NEGATIVE PHASE VELOCITY 



Oscillatory and excitable media are two kinds of media extensively investigated in non 



13j. CGLE is a typical 



linear dynamics. But antispirals are found only in oscillatory media 
and well known oscillatory system. Here we use ID CGLE with local pacing 17j| to explore 
nonlinear antiwaves. 
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We conduct simulations in an aray of N(N = 256) sites with Ax = 1.0 space discretization. 
Free boundary condition is used throughout the paper. Numerically, we add the periodic 
signal to the left boundary site of the ID system. We fix F — 1.0 for all simulations of this 
paper. The sign of the input frequency u in in Eq.(l) denotes the rotating direction of the 
forcing in complex plane, "+" clockwise and "-" anticlockwise. System (1) supports periodic 
travelling waves 



A = VT^¥e i ^ Uit+kx) (2a) 
u = a+{p- a)k 2 (2b) 

Equation (2b) together with the condition < k 2 < 1 define the constraint of u>, 

a<uj<fJ or f3<u<a (2c) 

In Fig. 1(a) we plot the input-output frequency relation[]j]] for (3 = —1.4, a = 0.2. There 
are four parameter domains in Fig. 1(a), corresponding to different response behaviors of 
the system to the local pacing. When u in < —0.41, the system cannot follow the pacing 
frequency due to too fast pacing frequency, the output frequency of the system is thus much 
lower than Ui n and also lower than the natural frequency (u = a), \uj ut\ < l^ol- We observe 
antiwaves asymptotically (see segment AB in Fig. 1(a) and the evolution pattern in t — x 
plane of Fig.l(b)). When —0.41 < uj in < 0, the system oscillates with the input frequency 
uj out = Uin resonantly and the waves propagate normally (CD segment in Fig. 1(a) and 
Fig. 1(c)). When < Ui n < 0.2, the system oscillates also with the input frequency u>i n while 
waves propagate inwardly to the wave source, i.e., one observes negative phase velocity (DE 




X X 



Fig.1 

FIG. 1: Numerical results of Eq.(l). a = 0.2, (3 = — 1.4, F = 1. (a) Input-output frequency 
relation. In AB and EF region, constraint Eq.(2b) forbids the complete drift 0J O ut = oj = oJi n , 
and the system motion is governed by the local dynamics, ui ou t ~ a and \k\ <C 1. CE region 
corresponds 1:1 resonant output io ou t = uji n . (b)-(f) Spatiotemporal patterns of the system for 
different parameter sets, (b) u>i n = —0.5 (AB region), antiwaves. (c) u>i n = —0.15 (CD region), 
normal waves, (d) u>i n = 0.15 (DE region), antiwaves. (e) LOi n = 0.21 (EF region). uj ou t ~ a = 0.2 
and k « 0, practically homogeneous oscillation, (f) uj out = u>i n = 0, stationary Turning pattern. 

segment in Fig.l(a) and Fig.l(d)). When uji n > 0.2, the system approaches asymptotically 
to nearly homogeneous oscillation, i.e. k 1, u out < o> = a (antiwaves with very large 
phase velocity, EF segment in Fig. 1(a) and Fig. 1(e)). An interesting point in Fig.l is that 
we can produce both NWs and AWs with the same parameter set of (a, (3) by changing the 
pacing frequency only. There is a special case of o>j„ = 0, which locates at the boundary 
between NWs and AWs. We expect that this local constant force produces neither normal 
waves nor antiwaves. Indeed we observe in Fig. 1(f) a pacing induced stationary Turning 
pattern. 

The above numerical results demonstrate that local pacing with suitable frequency can 
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produce AWs in nonlinear oscillatory systems indeed. These results show that nonlinear 
NWs and AWs are determined by the competition between the output frequency uj out and 
the local natural frequency u (u> = a for CGLE) of oscillatory media. In summary, if u out 
and ujq have the same sign while the absolute value of u out is smaller than that of u , we 
have 

u out ujQ > and|u;o| > \uj out \ Antiwaves (3a) 

In other two cases we have 

uJoutUjQ > and \u \ < \uj out \ Normal waves (3b) 
^out^o < Normal waves (3c) 

It is thus the competition between uj out and u that makes the system to select the correct 
sign of k in the dispersion relation Eq.(2b), and to determine positive (ku> > 0) or negative 
(kco < 0) phase velocity (or, say, normal or anti waves). 

It is emphasized that in Eq.(3) we use frequency relation between uj out and u , not u) in and 
ujq. Here oj out is the actual frequency of the system motion. In certain cases, the oscillatory 
system is completely driven by the pacing and we have oj out = uj in . Then uj in can be used 
instead of uj out for classifying the cases of Eqs.3(a), (b) and (c) (see regions CD and DE 
in Fig. 1(a)). In some other cases, the system cannot follow the rotation frequency of the 
pacing, and the asymptotic state of the system has frequency cu out ^ u in , then all the above 
analysis describes how the competition between u out and ujq produces normal and anti waves 
in oscillatory media. This is what happens in the segments AB and EF of Fig. 1(a). 



III. MECHANISM UNDERLYING NONLINEAR ANTIWAVES 

Now we try to understand the phenomena of negative phase velocity and antiwaves in 
CGLE. Firstly we begin with the analysis of the local dynamics of the system. Without the 
coupling and forcing terms the local dynamics of CGLE reads 

rl A 

— = A-(l + i a)\A\*A (4) 
Representing complex variable A by A = Re 10 , we can transform Eq.(4) to 



Eq.(5) is a typical limit cycle system with frequency a. Every grid of system (1) is a rotator 
e -iat w j 10 ro t a t; es clockwise (anticlockwise) in complex plane if a is positive (negative). 
Therefore, Eq.(l) is an oscillatory medium, and we call a the local natural frequency of the 
system, denoted as u = a. 

When we force the left boundary of the system, the grids near the left boundary are 
driven by the pacing, and they can quickly reach their asymptotic rotation with the output 
frequency u out . Due to the diffusion coupling, the perturbation from the left boundary can 
stimulate the grids distance away from the left boundary, and make these grids rotate with 
their natural frequency ujq in the first evolution stage. In the next evolution stage, the grids 
away from the pacing source will be synchronized to the asymptotic frequency u out , because 
the coupling restricts the phase difference between these grids and the grids near the left 
boundary. It is just the initial phase arrangement of different grids by u> ou t and uq plays the 
key role in determining the sign of the phase velocity. 

Now we analyze the three types of Eq.(3). First, we suppose that the pacing generated 
frequency u out is in the same direction with the local natural frequency ujq and we have 
\ujq\ > \u ou t\ (case Eq. (3a)). The grids near the left boundary rotate initially slower 
than its neighbor grids farther away from the pacing source. When synchronization of all 
these grids are reached the phase of the latter must be more in advance than that of the 
former. One thus observes phase propagation towards to the pacing source, an interesting 
phenomenon of negative phase velocity. This kind of phase setting is shown in Fig. 2 (a). 

Second, if u out and uo have the same sign while \u>\ > \u \ (case Eq.(3b)), the initial 
phase distribution of the grids is such prepared for the asymptotical travelling waves that 
any grid in the ID chain must have some phase delay in comparison with its neighbor 
grids nearer the pacing source as shown in Fig.2(b). In this case one can observe normal 
forward propagating waves. Third, if the two frequencies u out and u have opposite rotation 
directions (uJoUJ ut < 0, case Eq.(3c)), one can observes phase distribution of grids of Fig.2(c) 
yielding normal forward waves. 

Now we can give some clear conclusions on the problem of antiwaves and negative phase 
velocity. These conclusions are valid not only for CGLE systems, but also for general 
oscillatory systems. We suppose that a system is oscillatory with local dynamics having bulk 



frequency uj q 



13| and the medium supports periodic waves with frequency uj out generated 



by an arbitrary perturbation source. If uj out and ojq have the same sign and \uj out \ < \uq\, we 
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Fig. 2 



FIG. 2: Spatial phase distribution responding the perturbation from the wave source. Frequency 
uo is supposed in clockwise direction, (a) u out uJo > and \uj ou t\ > |^o|- Given a space site (with 
frequency to out), the sites farther away from the pacing, having frequencies (near loq) larger than 
uJouti rotate faster than the given site. When these sites are synchronized to the frequency io ou t, 
they must have phases more advanced than the phase of the given site. This phase distribution 
results in negative phase velocity and AWs. (b) uj out ujQ > 0, \co ou t\ > \ll>o\ an d (c) uj out uJ$ < 0. By 
the same argument as (a), the sites farther away from the pacing source have phase delay with 
respect to the given site. We can observe positive phase velocity and NWs. 

can surely produce antiwaves, and we can produce normal waves in all other cases, except 
two no generic critical cases: k = (producing homogeneous oscillations) and u out = 
(producing stationary Turning patterns). It is emphasized that our analysis based on the 
pictures of Fig. 2 does not use any special dynamic structure of CGLE, and thus does not 
need the requirement that the system is near Hopf bifurcation condition. This generality 
has never been explored in previous analysis on the problems of negative phase velocity. In 
our analysis, local oscillation (i.e., nonzero c^o) is a necessary condition for negative phase 
velocity. If the local dynamics of a medium does not have oscillation, we have ojq = 0. It is 
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clear from Eq.(3) that the system with ujq = can support only normally propagating waves. 
All excitable media have local frequency uq = 0, and we can therefore answer a previously 
raised question about antispirals in excitable media[3]: antiwaves (including antispirals) can 
never be found in excitable media. We have numerically computed CGLE for many other 
parameter sets, and computed different chemical reaction diffusion models with oscillatory 
and excitable local dynamics, the above conclusions have been fully confirmed. 

IV. PARAMETER DOMAIN FOR ANTIWAVES IN CGLE SYSTEM 

The analysis of Eq.(3) and the physical mechanism of Fig.2 are valid for general oscillatory 
systems. CGLE is a kind of oscillatory system universally appearing around Hopf bifurcation 
from a homogeneous stationary states to homogeneous oscillation. It is thus interesting to 
classify parameter domains for normal and anti waves of CGLE. In this section we focus on 
specifying the domain of antiwaves of CGLE in (a, j3, u ou t) parameter space. 

The motion of periodically forced CGLE can be periodic, quasiperiodic and chaotic. Here 
we are interested only in periodic travelling wave states. Far away from the pacing boundary, 
travelling waves must have solution of Eq. (2a) and the frequency u ou t must obey condition 
Eq.(2c). This condition is required for all periodic travelling waves, including both NW 
and AW. On the other hand, the conditions of Eq. (3) can be used to classify NW and AW. 
Eqs.(2c) and (3) together completely determine the parameter domains of NW and AW of 
CGLE. 

In Fig. 3(a) we present domain distribution of NWs and AWs in a— uo out plane with /3 fixed. 
In the parameter region marked "Forbidden" , the parameters do not satisfy the constraint 
of Eq.(2c), and no periodic travelling wave exists in this domain. In the blank region (the 
blow-left and up-right parts), the parameters satisfy Eq.(2c) and the conditions of Eq.(3a), 
we can observe AWs. In the shadowed region marked NW1 (NW2), the parameters satisfy 
Eq.(2c) and Eq.(3b) (and Eq.(3c)), and NWs are produced. In Figs. 3(b) and 3(c) we do 
exactly the same as Fig.3(a) with a and u fixed, respectively. Figs. 3(a), 3(b) and 3(c) are 
obtained for arbitrary positive (3, a and uj out . For negative (3, a and u out we can directly 
obtain domain distribution diagrams from Fig.3(a), 3(b) and 3(c), replacing a, u out ; (3, uj out ; 
and a, (3 by —a, —uj out ] —(3, —uj out ] and —a, —(3, respectively (e.g., see Fig.3(d)). 
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FIG. 3: Schematic figures of phase diagrams of CGLE, which are drawn, based on Eq.(2c) and (3). 
(a) Phase diagram in u} out — a plane with (3 > fixed. In the grey regions marked "Forbidden", 
travelling waves cannot be realized due to the condition of Eq.(2c). The blank region marked AW 
is antiwave region, characterized by Eq.(3a). The shadowed region marked NW1 (NW2) is a part 
of normal wave region, characterized by Eq.(3b) (by Eq.(3c)). (b) The same as (a) with a > 
fixed and the distribution of various characteristic domains is plotted in (3 — uw plane, (c) The 
same as (a) with uj out > fixed and the domain distribution is plotted in a — (3 plane, (d) The 
same as (c) with uj out < 0. (d) and (c) are symmetric against the symmetry center (a, (3) = (0,0). 

V. CONCLUSION 

In conclusion, we have found interesting phenomenon of negative phase velocity (anti- 
waves) induced by local pacing, and revealed the common mechanism underlying antiwaves: 
the competition between the actual system frequency uj out and the local natural frequency 
ujq selects the correct sign of the wave number k in the dispersion relation Eq.(2b) to pro- 
duce either normal waves or antiwaves. If condition uu > and \u\ < \u \ are satisfied, 
antiwaves can emerge definitely. To our best knowledge, it is the first time to emphasize 
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that both signs and values of to and u play essential roles in producing antiwaves. This 
understanding gives a convincing answer to the problem why no antispirals and antiwaves 
have been observed in excitable media. The criterion of antiwaves (Eq.(3a) and Fig. 2 (a)) is 
expected to be valid for general oscillatory systems, and this goes again beyond the previous 
results of negative phase velocity of CGLE. 
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